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ABSTRACT 

We report the discovery via radial velocity measurements of a short-period (P — 2.430420 ± 0.000006 
days) companion to the F-type main sequence star TYC 2930-00872-1. A long-term trend in the 
radial velocity data also suggests the presence of a tertiary stellar companion with P > 2000 days. 
High-resolution spectroscopy of the host star yields T c g = 6427 ± 33 K, \ogg = 4.52 ± 0.14, and 
[Fe/H]=— 0.04 ± 0.05. These parameters, combined with the broad-band spectral energy distribution 
and a parallax, allow us to infer a mass and radius of the host star of Mi = 1.21 ± 0.08 M Q and 
Ri = 1.09^Q'}g Rq. The minimum mass of the inner companion is below the hydrogen burning limit, 
however the true mass is likely to be substantially higher. We are able to exclude transits of the 
inner companion with high confidence. Further, the host star spectrum exhibits a clear signature of 
Ca H and K core emission indicating stellar activity, but a lack of photometric variability and small 
v sin / suggest the primary's spin axis is oriented in a pole-on configuration. The rotational period of 
the primary estimated through an activity-rotation relation matches the orbital period of the inner 
companion to within 1.5 er, suggesting that the primary and inner companion are tidally locked. If the 
inner companion's orbital angular momentum vector is aligned with the stellar spin axis as expected 
through tidal evolution, then it has a stellar mass of ~ 0.3 — 0.4 Mq. Direct imaging limits the 
existence of stellar companions to projected separations < 30 AU. No set of spectral lines and no 
significant flux contribution to the spectral energy distribution from either companion are detected, 
which places individual upper mass limits of M{ 2 ,3} ^ 1.0 M Q , provided they are not stellar remnants. 
If the tertiary is not a stellar remnant, then it likely has a mass of ~ 0.5 — 0.6 M Q , and its orbit is 
likely significantly inclined from that of the secondary, suggesting that the Kozai-Lidov mechanism 
may have driven the dynamical evolution of this system. 
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1. INTRODUCTION 

Exoplanet surveys have contributed to a wide range 
of ancillary astrophysical disciplines during the last two 
decades, including studies of variable stars, binary stars 
and brown dwarf (BD) companions. During the course 
of operation, these surveys detect a large variety of stel- 
lar binaries that can be used to study stellar structure, 
atmospheres and formation mechanisms. One example 
of the latter is a study of the multiplicity of close bina- 
ries, e.g., the fraction of close binaries that are in triple 
or higher-order systems. Indeed, triple systems are not 
uncommon amongst short-period binaries; 9 out of 16 bi- 
naries with P < 100 day s in the volume- limited sample of 
Ragh avan et al. I (|201Q[ ) are members of triple systems. 
Shorter-period binaries have an even greater probability 
of being in a multiple-star system (~ 80% for P < 7 days 
vs. - 40% for P > 7 days. iTokovinin et al. II2006T ). 

The orbital elements of such binaries, including the 
mutual inclinations of the companions' orbital angular 
momentum vectors, are fossil records of their formation 
process, and prov ide critical constraints to b inary star 
formation models (|Sterzik fc Tokovinin 1 12002). Compar- 
ison of the orbital and physical properties between dif- 
ferent binary hierarchi es also provides in sight into binary 
star formation theory ()Tokov inin 2008). In fact, the dy- 
namical evolution of these systems may be dominated by 
dynamical interactions between the inner and outer com- 
panions via a combination of the Kozai-Lidov mechanism 
(|Kozai Ifl962t iLidov II1962D and tidal forces, which drive 
the inner companion to shorter orbital separations until 
it circularizes with some period P < 10 days, beyond 
which tidal forces are ineffective (Fabrvc kv fc Tremaine I 
l200l . 

In this paper, we present the discovery of a compan- 
ion with a substellar minimum mass orbiting the bright 
(V = 9.8) F-type star TYC 2930-00872-1 faog et al. I 
120001 hereafter TYC 2930), with an orbital period of 
P = 2.430420 ± 0.000006 days. This discovery is part 
of a series of papers dedicated to analyses of individual 
low-mass companions in anticipation of a global analysis 
of the MARVELS sample at the conclusion of th e survey 
(e.g.. iLee et al. 1 120111: [Wisniewski et al. 1 120121 ), there- 
fore, TYC 2930 is also designated "MARVELS-2" as an 
internal reference within this series. The a priori transit 
probability of the inner companion is ~13% with an ex- 
pected central transit depth of ^0.9 ± 0.25% for a 1 Rj up 
companion radius, although no transits are detected. An 
additional, long-term trend in the RV data is detected 
from a stellar tertiary in the system. A detailed analysis 
of the combined radial velocity, spectroscopic and pho- 
tometric data suggests the inner companion is oriented 
towards a pole-on configuration, and is more likely an M 
dwarf with a mass ^0.3 — 0.4 M Q , while the tertiary is 
likely to be less inclined. In such a scenario, the mutual 
inclination between the secondary and tertiary is likely 
to be significant, which would make this an excellent ex- 
ample of a system whose dynamical history was driven 
via the Kozai-Lidov mechanism. 

The paper is organized such that £12.11 describes the 
spec troscopic observations and their data processing, 
£12.21 describes the archival and observed photometry for 
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the system, describes the characterization of the host 
star's properties, including mass, radius, effective tem- 
perature, surface gravity, metallicity, stellar activity and 
rotation rate, §H describes our determination of the or- 
bital parameters from fitting the measured RVs, £|5] de- 
scribes both Lucky Imaging and adaptive optics imaging 
to search for any wide companions to TYC 2930, fJ5] de- 
scribes our search for photometric variability and any po- 
tential transits of the inner companion, £J7] discusses the 
tidal evolution of the inner companion, £J5] describes the 
posterior distribution of the true masses for both the sec- 
ondary and tertiary given the results from the previous 
sections, and finally, £J5] investigates the possible dynami- 
cal history of the system via the Kozai-Lidov mechanism. 

2. DESCRIPTION OF OBSERVATIONS 
2.1. Spectroscopic Observations 

MARVELS (Multi-object APO Radial Ve l ocity Ex- 
oplanets Large-area Survey, iGe et al.~l I2008D is 
part of the Sloan D i gital Sky Survey III (SDSS- 
III, Eiscnste in et al. I 1201 ID . The instru ment uses 
dispe r sed fixed-d e lay inter f erome t ry (DFDI: IGe et al. I 
2003 [Gel [2001 lErskine I [200l lErskine et al. I 120031: 
van Ev kcn et aTTI 1201 ID on the 2.5m SDSS telescope 
(jGunn et al .1120061 ) at Apache Point Observatory (APO) 
to measure precision radial velocities of 60 stars simul- 
taneously. Both beams of the interferometer are imaged 
onto the detector for a given star, for a total of 120 spec- 
tra, producing two simultaneous RV measurements for 
each star from beams that travel through a slightly dif- 
ferent instrument path. The survey began in the Fall of 
2008 and will ultimately target several thousand stars be- 
tween 7.6 < V < 12, with a baseline goal of < 30 m s _1 
precision for the faintest stars. Each star is observed 
^20-30 times over a typical baseline of 1.75 years. In ad- 
dition to exoplanets, the survey will conduct studies of 
stellar atmospheres, binary stars, and rare companions 
such as BDs and very low mass (VLM; M < 150 Mj up ) 
stars at short orbital periods. 

TYC 2930 was observed a total of 33 times over a base- 
line of 707 days. The data were processed by the MAR- 
VELS pipeline following the steps described in lLee et al. I 
(|2011l ). The resultant RV measurements from both inter- 
ferometer output beams were combined via a weighted 
average after they were found to agree to within the 
measurement errors. The formal mean RV precision was 
23 m s" 1 . Following iFleming et al~l (|2010f l . the RV un- 
certainties for this star were further scaled by a "quality 
factor" QF = 6.69, a first-order correction used to par- 
tially account for residual systematic errors. For each of 
the other 118 spectra in this field, an individual QF is 
calculated as the RV root-mean-square (rms) about the 
mean, divided by the median formal RV uncertainty for 
that star. On average, most of the MARVELS targets 
in a given field should be RV-stable at the level of tens 
of m s _1 , and therefore should have QF ~ 1. Since the 
average QF across the plate is significantly larger, we 
treat that as one measurement of the residual uncertain- 
ties from the pipeline-produced RVs. The dates and RVs 
from the MARVELS observations are presented in Table 

m 

Additional RV observations were conducted using the 
Spettrografo Alta Risoluzione Galileo (SARG) spectro- 
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graph (|Gratton et al. 1 120011 ) on the 3.58m Telescopio 
Nazionale Galileo (TNG) telescope. The data were ob- 
tained using the yellow grism with a slit of 0.8" x 5.3" on- 
sky that produces a resolving power of R — 57000 over 
the wavelength range 462 < A < 792 am. A total of 20 
observations were taken, spanning ~408 days, using an 
iodine cell that serves as an RV calibration source. The 
average signal-to-noise ratio per resolution element, aver- 
aged across the central 200 pixels of all the orders, ranges 
from 150-290. An additional observation was taken with- 
out the iodine cell to be used as a stellar template and to 
derive stellar parameters. The signal-to-noise ratio per 
resolution element of the template spectrum is ~400 at 
607 nm. 

The SARG data are processed using the standard 
IRAF Echelle reduction packages. Frames are trimmed, 
bias subtracted, flat-field corrected, aperture-traced and 
extracted. ThAr lines are used to calibrate the wave- 
length solution. The RVs are measured using the iodine 
cell technique (jMarcv fc Butler I H99l . The 21 SARG 
orders that have sufficiently strong iodine lines lie in the 
wavelength range 504 < A < 611 nm. Each order is sub- 
divided into 10 sections from which an RV is measured. 
The resulting 210 RV measurements are then 2-cr clipped 
using three iterations. The remaining N RV measure- 
ments are averaged to produce a single RV measurement. 
The RV uncertainty is given by ctrv = °~ -/V -1 / 2 , where 
cr is the dispersion of the points after the 2-tr clipping. 
Table [2] contains the dates and RVs from the SARG ob- 
servations. 

A high signal-to-noise ratio spectrum of TYC 2930 was 
obtained with the ARC Echelle Spectrograph (ARCES; 
iWang et al~ll2003D on the APO 3.5m telescope for the 
purposes of stellar characterization. The spectrograph 
delivers R ~31500 spectra spanning a wavelength range 
320 < A < 1000 nm on a single 2048x2048 SITe CCD. 
The spectra were reduced using an IRAF script that cor- 
rects for bias and dark current subtraction, cosmic rays 
and bad pixels. Flat-Fielding is performed using a com- 
bination of quartz lamp exposures with and without a 
blue filter, while a ThAr lamp is used for wavelength cal- 
ibration. A single integration of 900 seconds was taken, 
yielding a spectrum with a signal-to-noise ratio of ^220 
per resolution element at 607 nm, 

Long-term, queue-scheduled iSliot rone et al. I I2007D . 
RV monitoring of the TYC 2930 system has been ini- 
tiated using the High Resolution Spectrograph (HRS, 
Tulll 119981) on the Hobby-Eberly Telescope (HET, 
Ramsey et al. IIT998I ) to further characterize the orbit of 
the suspected long-period companion. These RVs are 
expected to be presented in a separate paper at the con- 
clusion of that project. 

2.2. Photometry Observations 

Photometry of TYC 2930 was performed using the 
Hereford Arizona Observatory (HAO), a private ob- 
servatory in southern Arizona (observatory code G95 
in the IAU Minor Planet Center). Observations 
were taken in Johnson B and V filters using a 
Meade 14-inch LX200GPS telescope and a 2184x1472 
pixel SBIG ST-10XME CCD. Landolt standard stars 
(jLandolt fc Uomoto 1120071 lLandolt 1 12009) were observed 
in the Kapteyn Selected Area 98 (SA 98) for calibration. 
A photometric precision of 0.023 mag was obtained in B 



and 0.018 in V. The measured fluxes and uncertainties 
are presented in Table [3j 

We obtained relative photometric time series from 
several ground-based telescopes (SuperWASP, Allegheny 
Observatory, KELT-North) to search for transits and ex- 
amine the photometric stability of the primary star. We 
briefly describe each of these data sets in turn. The Su- 
perWASP instruments measure fluxes of millions of stars 
via wide-angle images of the night sky using a broad- 
band filter that covers 400-700 nm, and are described in 
iPollacco et al~l ()2006l ). For TYC 2930, a total of 2204 
observations from 2006 and 1309 observations from 2007 
are extracted from the SuperWASP public archivj 2 ^!. 

We obtained photometric observations on 7 nights 
in February and March of 2011 using the Keeler 16- 
inch Meade RCS-400 telescope at Allegheny Observatory. 
The CCD detector is a 3060x2040 pixel SBIG KAF- 
6303/LE with a 0.57" per pixel scale, and all observations 
were taken through a Johnson-Cousins R filter. Typi- 
cal seeing was 2.5" with integration times ranging from 
20-30 seconds. The images were processed using stan- 
dard bias, dark, and flat-field calibration images taken on 
the same nights. Astrometric solutions were computed 
based on the positions of stars in the 20' x 30' field of view 
from the 2MASS Point-Source Catalog (jSkrutskie et al. I 
2006). After image calibration, we performed circular 
aperture photometry with a 10-pixel radius (5.7" on sky) 
and estimated the local sky background from a 15-20 
pixel annulus around each star. Relative photometry was 
determined by comparing the measured flux from TYC 
2930 with two nearby stars. 

We also extracted photometric time series data of TYC 
2930 obtained by the Kilodegree E xtremely Little Tele- 
scope (KELT) Nort h transit survey (jPepper et al. 1120071 : 
ISiverd et al. 1 12009?) . KELT uses a red-pass filter with a 
50% transmission point at 490 nm, which, when folded 
with the CCD response, yields an effective bandpass sim- 
ilar to R, but broader. The standard KELT data reduc- 
tion procedure uses the ISIS image subtraction package 
(jAlard fc Lupton"1 [l998l) , combined with point-spre ad fit- 
ting photometry using DAOPHOT (jStetson II1987T ). 

In the case of TYC 2930, the standard KELT data re- 
duction procedure yielded an unusable lightcurve due to 
the presence of the nearby bright star HD 42903, which 
was partially saturated in the KELT images. We correct 
the systematics by performing simple aperture photom- 
etry on both TYC 2930 and HD 42903 using the sub- 
tracted images. We used two apertures centered on HD 
42903, and one aperture centered on TYC 2930. We sized 
the apertures around HD 42903 such that they formed 
an annulus that included the systematic artifacts. The 
single aperture around TYC 2930, in the middle of the 
artifact, had the same diameter as the width of the an- 
nulus around HD 42903. By subtracting the summed 
flux in the aperture around TYC 2930 from the annulus 
around HD 42903, we are left with the average negative 
flux value in the artifact for each of the subtracted im- 
ages. We then used this average value to correct the 
results from the aperture around TYC 2930. This proce- 
dure was tested using known variable stars and on stars 
with similar brightness ratios and angular separations to 
confirm that accurate results were obtained and no in- 

23 http://www.wasp . le . ac.uk/public/ 
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trinsic variations were suppressed. 

3. STELLAR CHARACTERIZATION 
3.1. Stellar Parameters 

TYC 2930 (HIP 29714) is a bright F-type star located 
6.16' from the center of the open cl uster NGC 2192. The 
Hipparcos parallax measurement (jvan Leeuwen 1 120071 ) 
places t he star at a distance of d = 1 40 ± 29.5 pc. The 
RPM-J (jCollier Cameron et al. Il2007| ) value of 2.14 and 
(J — H) color of 0.23 are consistent with a main sequence 
star. We further characterize the host star's proper- 
ties using the SARG template spectrum and the ARCES 
spectrum by measuring equivalent widths (EWs) of Fe I 
and Fe II lines. We utilize two independent pipelines 
that derive stellar atmospheric parameters based on the 
Fe I and Fe II excitation and ionization equilibria. We 
refer to these different pipelines as the "IAC" (Instituto 
de Astrofisica de Canarias) and "BPG" (Brazilian Par- 
tic ipation Group) pipelines, which are described in detail 
by IWisniewski et al. I ((2012). We apply both these anal- 
yses to the SARG and ARCES spectra, and find spectro- 
scopic parameters that are consistent across both groups 
and both instruments. We summarize the individual T c g , 
log (g), [Fe/H], and & in Table H 

A fin al, mean value for each pa rameter is calculated fol- 
lowing Iwisniewskreran ()2012l ). yielding T cB = 6427 ± 
33 K, log (5) = 4.52 ± 0.14, [Fe/H] = -0.04 ± 0.05 and 
£t = 1.40±0.05 km s _1 . We note that while there can be 
correlations between the measured stellar parameters, we 
treat their uncertainties as independent and Gaussian- 
distributed in this analysis. Estimates of the primary's 
mass and radius are determined using a Markov chain 
Monte Carlo ( MCMC) analy s is app lied to the empirical 
relationship of iTorres et al. I (|2010l ). the Hipparcos par- 
allax, and the stellar parameters described above. The 
uncertainties for the mass and radiu s include the correla - 
tions of the best-fit coefficients from lTorres et al~1 (j2010f ) 
and the reported scatter in that relation (<7i ogTO ~ 0.027 
and CTi g r = 0.014). The radius of the primary is 



+0.15 
-0.13 ^0 



Rn and the mass is M = 1.21 ± 0.08 M 



0- 



R = 1.09 

All of the stellar parameters are summarized in Table [H 

3.2. Rotation Rate, SED Fitting, and Stellar Activity 

We measure the stellar rotational velocity v sin / from 
the SARG template spectrum. Note that we utilize a no- 
tation that distinguishes / (the angle between our line- 
of-sight and the stellar rotation axis) from i (the angle 
between our line-of-sight and a companion's orbital angu- 
lar momentum ve ctor). We use an interpolated Kurucz 
model spectrum (|Kurucz Ill993|) using the spectroscop- 
ically determined T c g, log (g) and [Fe/H], convolved to 
the instrumental profile (FWHM of 5.3 km s _1 ). Testing 
showed that macroturbulence (Q) had only a marginal 
affect on the final result, so we adopt values ranging 
from 2-5 km s -1 . The model spectra are broadened 
with a Gray's profile over a range of v sin / from 0- 
10 km s . These models are then compared via a 
X 2 analysis with the observed spectrum, which yields 
usin/ = 3.8 (+1.9, —2.8) km s _1 , effectively placing an 
upper limit of vsinl < 6 km s . 

We construct an SED using fluxes from GALEX 
(|Martin et al. II2005D. t he HAO observations, the 2MASS 
Skrutskie et al. I [21)01 ii Point Source Catalog, and the 



four WISE (IWright et al. l[2010l) bands. NextGen models 
from lHauschildt et al.l ( 1999f ) are used to construct theo- 
retical SEDs by fixing T c s , log (g) and [Fe/H] at the spec- 
troscopic values, while the extinction Ay is constrained 
to a max imum value of Ay = .6 based on the reddening 
maps of iSchlegel et al.l (|1998f ) for galactic coordinates 
(/,&) = (173.365948°, 10.729936°). Fig. □ shows the best- 
fit model, which has a x 2 / dor = 1.2, A v = 0.33±0.06, no 
evidence for IR excess, and some excess in the GALEX 
FUV band. Fig. [5] places the star o n an HR diagram 
based on Yonsei-Yale stellar models (jDemarque et al.~l 
I2004D . indicating that TYC 2930 is consistent with an 
F-type dwarf with an age t < 2 Gyr. 

To further explore the FUV excess, Fig.[3]compares the 
ARCES spectrum of TYC 2930, centered on the Ca II 
K line at 393.37 nm, with archival Fibre -fed Extended 
Rang e Optical Spectrograph (FEROS, iKaufer et al. I 
119991) spectra of the s tandard stars HP 43 042, HD 142 
and HD 120136 from IGhezzi et al. I (pOlOl ). FEROS is 
an R ~ 48000 spectrograph with a wavelength range of 
350 < A < 920 nm and high throughput (~ 20%). There 
is clear Ca K core emission from TYC 2930 indicating 
signifi cant chromospheric activity. We measure the S 
index (IVaughan et al. Ill978t IVaughan fe Preston Ill980t 
iDuncan et al. 11199 If ) from the APO spectrum and con- 
vert to R' HK , finding a \ogR' HK = —AAA ± 0.05. 

4. ORBITAL ANALYSIS AND COMPANION MINIMUM 
MASS 

4.1. Radial Velocity Fitting 

The MARVELS RVs show evidence of a long-term, 
positive linear trend indicating a possible tertiary ob- 
ject. The TNG data also showed evidence of a long- 
term trend, but with a negative slope. Fitting the MAR- 
VELS+SARG data with a single-companion model com- 
bined with non-Keplerian trends (linear, parabolic, cu- 
bic) yielded residuals with significant systematics, sug- 
gesting a two-companion Keplerian model is required. 

The co mbined RVs were initial ly fit using the RVLIN 
package ( Wri ght fe How ard 2009]) for the purpose of ob- 
taining initial values of the orbital parameters. Un- 
certainties are calculated later using MCMC analysis. 
The initial best-fit orbital period for the inner com- 
panion is P2 = 2.430420 days, with a semi-amplitude 
K 2 = 8724 m s _1 and an eccentricity that is consistent 
with a circular orbit. Fig. 0] shows the MARVELS (blue) 
and SARG (red) RVs phase-folded on the best-fit orbital 
solution for this inner companion after removing the ef- 
fects of the longer period orbit. The bottom panel plots 
the residual RVs after removing the shorter period orbit. 
The unfolded RVs from MARVELS (blue) and SARG 
(red) are shown in Fig. [5j where the shorter period orbit 
has been removed and the best-fit model of the longer 
period orbit is plotted as the solid line. The residuals of 
the combined, two-companion solution are shown in the 
bottom panel. 

To derive final orbital parameters and associated un- 
certainties, we perform an MCMC a nalysis closely fol- 
lowing the methods of iFord I ((2006) . For review, our 
goal is to estimate the uncertainties in our set of model 
parameters, 8 = {P 2 , P 3 , K 2 , K 3 , e 2 , e 3 , u> 2 , w 3) T P2 , 
Tpz, 7off, 7o,inst}, where P is the orbital period, K is 
the RV semiamplitude, e is the orbital eccentricity, oj 
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is the argument of periastron, T p is the epoch of peri- 
astron, 7 g is the offset between the two sets of instru- 
ments, 7o,inst is the (instrumental) systemic velocity, and 
the subscripts j = {2, 3} refer to the shorter period and 
longer period companions, respectively. We sample the 
posterior probability distribution given by B ayes' theo- 
rem, where specifics on the priors and likelihood fun ction 
can be found in Section 3 of Zakamsk a et al. I (|201 If ) . To 
help accelerate convergence, we use additi onal combina - 
tions of parameters identified in Section 4 of lFord I {2006). 
We do not attempt to place constraints on stellar jitter 
in our model. 

We test for non-co nvergence by monit oring the 
Gelman- Rubin statistic (|Gelman et al. I [20031 ). verifying 
it is less than 1.02 for each of the parameters, and that 
chains have been allowed to run long enough to enter 
these regions of parameter space at least 100 times. The 
orbital parameters and l-er-equivalent confidence levels 
for the inner companion are given in Table [5j The outer 
companion's orbital parameters are not well-determined, 
since the orbital period is longer than the baseline of the 
measurements, but we place a lower limit of P3 > 2000 
days. The RV semiamplitude and eccentricity of the 
outer companion are positively correlated with the best- 
fit orbital period, which must be accounted for when con- 
straining the outer companion's properties. 

To place the RVs on an absolute scale, SARG spectra 
from 620-800 nm are cross-correlated with a high reso- 
lution solar spectrum. To partially account for temporal 
variation in the slit illumination and wavelength solution, 
a correction is applied via cross-correlation of the telluric 
lines near 690 nm with a nu merical mask. The telluri c 
line locations are taken from iGriffin fc Griffinl (fl~973l ); 
the corrections are typically a few hundred m s _1 . Af- 
ter removing the barycentric velocity and orbital motion 
of the companions based on the MCMC parameters, we 
find a median absolute RV of 70 = 35.751±0.285 km s _1 , 
where the uncertainty is taken as the rms about the me- 
dian. 

The {U, V, W} space velocities can then be calcu- 
lated using the absolute RV along with the paral- 
lax and proper motion measurements from Hippar- 
cos (Ivan Leeuwen I [2001 . We find {U,V,W} = 
{-31.2 ± 0.7, -7.4 ± 3.1, 3.0 ± 1.6} km s" 1 , where U is 
pointing towards t he Galactic center. F rom the classi- 
fication scheme of IBensbv etaTl ((2003T ). TYC 2930 is 
almost certainly a member of the thin disk, as its rel- 
ative probability of being a thick disk member is just 
0.71 ±0.02%. 

4.2. Mass Functions of the Secondary and Tertiary 

Using the MCMC chain from the joint RV fit, we can 
derive the mass functions Mj of companion j = 2, 3, 



= (M 3 sing 3 
3 ~ (Mi + Mj) 2 



2nG 



(1) 



The mass functions are the only properties of the com- 
panions that we can derive that are independent of the 
properties of the primary. For the secondary, we find, 



M 2 = (1.6711 ± 0.0050) x 10~ 4 M G 



(2) 



where the uncertainty is essentially dominated by the 
uncertainty in K 2 , such that om 2 /^2 ~ ^uk^/ K%) = 



3 x 0.1% - 0.3%. 

For the tertiary, the uncertainty in the mass function is 
much larger, because of the incomplete phase coverage of 
the radial velocity curve (Figs. [4] and [5]). In particular, 
there is a broad tail toward high mass functions. We 
therefore quote the median and 68% confidence interval, 



M 3 = 2.55tlH x 10~ 2 M o . 
4.3. Minimum Mass and Mass Ratio 



(3) 



To determine the mass or mass ratio of the secondary 
and tertiary, we must estimate the mass of the primary, 
as well as the inclination of the secondary and tertiary. 
To estimate the mass and radius of the primary, we use 
an MCMC chain where, for each link in the MCMC chain 
from the joint RV fit, we draw a value of T e g, log g, and 
[Fe/H] for the primary from Gaussian distributions, with 
me ans and dispersions given in Table [3] We then use 
the iTorres et al. I (|2010D relations to estimate the mass 
Mi and radius R± of the primary, including the intrinsic 
scatter in these relations. 

The minimum mass (i.e., Mi if sin 12 = 1) and mini- 
mum mass ratio of the secondary are: 

M 2 , min = 68.1 ± 3.0 M Jup = 0.0650 ± 0.0029 M Q , (4) 

q 2 = 0.0535 ±0.0012 (5) 

The uncertainties in these estimates are almost entirely 
explained by the uncertainties in the mass of the primary: 
o-m 2 /M 2 ~ (2/3)(<WMi) = (2/3) x 6.7% ~ 4.5%, al- 
most exactly the uncertainty in M 2 m i n above (4.4%), 
and ajq ~ (l/3)(a Ml /Mi) ~ (1/3) x 6.7% ~ 2.3%, the 
uncertainty in q (2.3%). As we show in $6.31 an edge- 
on orbit for the secondary is excluded from the lack of 
transits for reasonable assumptions about its radius. 

The minimum mass and mass ratio of the tertiary are 
much more poorly constrained due to the incomplete 
phase coverage of the orbit. We find median and 68% 
confidence intervals of, 



M 



3, min 



426 



+261 
-98 



M 



Jup 



0.407™ M , (6) 
(7) 



q 3 = 0.334] 



0.0761 



The uncertainties in these quantities contain significant 
contributions from both the uncertainty in the host star 
mass and the tertiary mass function. 

5. IMAGING 



Lucky Imaging (Fried] I1978D was performed in Oct . 
2010 and Oct. 2011 using FastCam (|Oscoz et al. 1 120081) 
on the 1.5 m TCS telescope at Observatorio del Teide 
in Spain to search for companions at large separations 
from the primary star. Lucky Imaging consists of tak- 
ing observations at very high cadence to achieve nearly- 
diffraction-limited images from a subsample of the total. 
During the Oct. 2011 observations, the CCD gain was 
adjusted, and therefore that night's data are analyzed as 
two different image sets. For the Oct. 2010 data, a total 
of 140000 frames comprised of 50 ms integrations were 
obtained in the I band spanning 21"x21" on sky. For 
the Oct. 2011 observations, a total of 31000 frames com- 
prised of 50 ms integrations were obtained in the low-gain 
setting, and 100000 frames comprised of 40 ms integra- 
tions were obtained in the higher-gain setting. Image se- 
lection is applied using a variety of selection thresholds 
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(best X%) based on the brightest pixel (BP) method, 
making sure that non-speckle features are avoided. 

The BPs of each frame are then sorted from brightest 
to faintest, and the best X% are then shifted and added 
to generate a final image, where X = {15, 85}. The effec- 
tive Strehl ratios for X = 85% are {0.036, 0.037, 0.043} 
for the three image sets, respectively. Fig. [5] shows com- 
posite images for X = 15 and 85% of the frames for each 
set of observations. The intensities are detector counts 
on a linear scale after being stacked and normalized by 
the number of images used in the stacking. The artifact 
in the Oct. 2010 frames is a result of imperfect telescope 
tracking. No companion is detected a t the 3-cr level, 
where a is defined using the procedure in lFemema et al.1 
poll based on the rms of the counts within concentric 
annuli centered on TYC 2930 and using 8-pixel boxes. 

In addition to the Lucky Imaging, we conducted adap- 
tive optics (AO) imaging to search for any wide stellar 
companions to TYC 2930. The Keck AO images were ob- 
tained with NIRC2 (PI: K. Matthews) on 2011 August 
30 UT. The observations were conducted in the K' band 
using the narrow camera setting, resulting in a plate scale 
of 9.963 mas pix" 1 (|Ghez et al. 1 12008( 1 . The total inte- 
gration time was 65 sec using a three-point dither pat- 
tern. Images were processed using standard pixel clean- 
ing, flat-fielding and stacking procedures. Fig. [7] shows 
the processed Keck AO image; no evidence for wide 
stellar companions can be seen. Detectability curves 
(3-cr) are calculated as a function of separation from 
TYC 2930 for both the Lucky Imaging and AO data. 
Contra st ratios are converted into mass sensitivities us- 
ing the IBaraffe et al.~l (I2003D models for the Keck band 
and IGirardi et al. I ( 20021) models for the Lucky Imag- 
ing band. As can be seen in Fig. [5J we can exclude 
stellar companions at projected separations greater than 
^50 AU. While the Keck constraints are superior com- 
pared to the Lucky Imaging constraints, they also rely on 
a very expensive resource (namely, the Keck telescope). 
Since Lucky Imaging can be conducted on much smaller 
and more readily available telescopes, it is a good re- 
source to use in the search for wide companions in the 
absence of 10-meter telescope access. 

6. ANALYSIS OF THE RELATIVE PHOTOMETRY 
6.1. Summary of Datasets 

The WASP photometric dataset for TYC 2930 con- 
sists of 3975 points spanning roughly two years from 
HJD' = 3831 to 4571. The full, detrended WASP dataset 
has a relatively high weighted rms of 2.9% and exhibits 
evidence for systematics. The distribution of residuals 
from the weighted mean is asymmetric and highly non- 
Gaussian, showing long tails containing a much larger 
number of > 3cr outliers than would be expected for a 
normally-distributed population. 

We clean the WASP data by adding in quadrature to 
the photon noise a systematic uncertainty (cr sys ) that re- 
sults in a distribution of residuals closest to the Gaussian 
expectation. We reject the largest, error-normalized out- 
lier from the mean flux value and scale the uncertainties 
by a factor r to force x 2 /dof = 1, iterating until no more 
outliers > Aa remain. Although 4cr is a slightly larger de- 
viation than we would expect based on the final number 
of points, we adopt this conservative threshold to avoid 



removing a potential transit signal. We find r — 0.39 and 
(T S y S = 0.0053, retaining 3731 data points with an rms of 
0.71% and x 2 /dof = 1 (by design). 

The Allegheny photometric dataset consists of 1280 
points spanning roughly 44 nights from HJD'=5596 to 
5640. The weighted rms of the raw light curve is 0.48%; 
this is a factor of ~ 3 times smaller than the average frac- 
tional photometric uncertainties, indicating that these 
errors have been overestimated. Although there is no 
clear evidence for systematic errors in this dataset, we 
repeat the identical procedure as with the WASP data 
for consistency. We find r = 0.38 and a sys = 0.0011, 
with a final rms of 0.42% from 1274 data points. 

The KELT dataset consists of 2781 data points span- 
ning roughly 3 years from HJD' = 4107 to 5213. The 
weighted rms of the raw light curve is 0.62%, and the 
mean uncertainty is 0.55%, indicating that these are rea- 
sonably well estimated. Nevertheless, for consistency we 
clean the data in the same way as the other two data sets, 
finding no outliers > 4cr, r = 1.09, and er sys = 0.0021, 
with a final rms of 0.62%. 

Finally, we combine all the relative photometry af- 
ter normalizing each individual data set by its mean 
weighted flux. The top panel of Fig. [9] shows the com- 
bined data set, which consists of 7786 data points span- 
ning roughly 4.4 years from HJD' = 4022 to 5640, and 
has a weighted rms of 0.58%. The resulting light curve 
is constant to within the uncertainties over the entire 
time span. Within the KELT data set, which spans ~ 3 
years, we find no strong evidence for long-term intrinsic 
variability at a level > 0.6%. 

6.2. Search for Periodic Variability 

We ran a Lomb-Scargle periodogram on the full 
dataset, testing periods between 1 — 10 4 days. The re- 
sulting periodogram, shown in the bottom panel of Fig. 
[51 displays a large number of formally significant peaks. 
The inferred amplitudes are all < 0.1%; similar to the 
level of systematic errors inferred when cleaning the light 
curves. A comparison of periodograms performed on 
the individual data sets demonstrates that the strongest 
peaks arise from only one dataset, and are not corrobo- 
rated by the other datasets. The most significant peak 
in the combined dataset has a period of 8.24 days, with a 
power of ~ 100 and an amplitude of ~ 0.13%. However, 
as shown in the inset, the signal comes almost entirely 
from the WASP dataset (blue). In general, the KELT 
dataset (red) shows significantly reduced power on all 
periods < 100 days; the rms of the periodogram in this 
range is only ~4, as compared ~14 for the combined data 
set. Although the different results inferred for different 
datasets could in principle arise from real variability that 
is not strictly periodic or persistent, it is more likely that 
there exist systematics in the data sets in the form of 
residual correlations on a range of timescales. 

Restricting attention to periods within 10ct of that in- 
ferred for the companion (P = 2.430420±0.000006days), 
the maximum power is ^20 with an amplitude of only 
^0.06%, a factor of ~10 times smaller than the rms of 
the light curve. Considering an expanded range of pe- 
riods within 2<T of the period of the primary as inferred 
from the R' HK index (P = 2.93 ± 0.37 days, see $7$, the 
maximum power is ~46 with an amplitude of ~ 0.09%. 
We do not regard these maxima as significant, and con- 
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elude that the star does not exhibit periodic variability 
at either the expected rotation period of the star or the 
period of the inner companion at a level > 0.1%. 

Fig. [10] shows the combined light curve, folded at the 
median period and time of conjunction of the companion, 
(P = 2.430420 ±0.000006 days and T c = 2454842.2640), 
and binned in phase using bins of 0.05. The weighted 
rms of the binned light curve is ~ 0.087%. Although 
the variations are larger than expected from a constant 
light curve based on the uncertainties (x 2 /dof ~9), we 
again suggest that these are due to systematic errors in 
the relative photometry. In particular, the folded, binned 
KELT light curve (red) shows a somewhat lower rms of 
~ 0.068% with a x 2 /dof = 1.8, and is more consistent 
with a constant flux. We conclude that there is no strong 
evidence for variability of TYC 2930 on any time scale we 
probe. We can robustly constrain the amplitude of any 
persistent, periodic variability to be less than < 0.1%, 
and we can constrain the amplitude of photometric vari- 
ability at the period of the companion to < 0.07%. 

Given the estimated stellar mass, the companion pe- 
riod, and the minimum mass, the amplitude of ellipsoidal 
variability is expected to be: 

S cmp ~ 0.03% (sin/) ( m2 SinI ) ( -J*J) ( -?-) ' 
P V ' \67 M Ju J ^1.2 Mq J v 2.43 dj 

(8) 

with a period of P/2 (jPfahl et al. II2008D . This expected 
signal is compared to the binned data in Fig. 1 1 01 demon- 
strating that it is just below the level of detectability. 
Since smaller inclinations lead to lower amplitudes, we 
are unable to constrain the inclination using ellipsoidal 
viability. 

6.3. Excluding Transits of the Secondary 

The probability that a low-mass companion transits 
can be determi ned given the orbital par ameters from 
the RV solution (|Kane fc von Braun II2008T ). Assuming a 
uniform distribution in cos i, the a priori transit probabil- 
ity for the secondary is relatively high, ~13%. However, 
the light curve folded on the ephemeris of the inner com- 
panion shows no evidence for a transit at the expected 
time of conjunction, with an upper limit to the depth of 
any putative transit of < 0.2%. In contrast, the central 
transit of a Jupiter-sized companion would be expected 
to have a depth of S ~ (r/Ri) 2 ~ 0.9% and a duration of 
~ 0.042 in phase. We conclude that our observations rule 
out such a transiting companion with high confidence. 

We perform a quantitative search for tr ansit signals us- 
ing a m ethod similar to that described in lFleming et aLl 
(2010). We use the distributions of the secondary period 
P, semiamplitude K, and time of conjunction T c from 
the MCMC RV analysis, setting the eccentricity of the 
secondary to zero for simplicity. For each link in the 
MCMC chain, we draw a value of T e g, log g, and [Fe/H] 
for the primary from Gaussian distributions, with means 
and dispersions given in Table [3] and dete rmine the pri- 
mary' s mass Mi and radius R\ using the iTorres et al. I 
(2010) relation on these values. 

We then draw a value of cos i from a uniform dis- 
tributiorF^I. and use the resulting values of P, K, Mi, 

24 Formally, this assumes a prior on the companion mass m p 
that is uniform in logm p . 



J?i and i to determine the secondary mass m p , semima- 
jor axis a, and impact parameter of the secondary orbit 
b = acosi R^ 1 . Finally, adopting a radius for the com- 
panion, r, we determine if the companion transits, and 
if so we determine the properties of the l ight curve us- 
ing the routines of [Mandel fc Agofl (|2002T ) . We assume 
quadratic limb darkening and adopt coefficients appro- 
priate for the R band from lClaret fc Hauschildt I (|2003|L 
assuming solar metallicity and the values of T e g and log g 
listed in Table [3] For reference, Fig. [10] shows the pre- 
dicted transit signatures for the median values of the 
physical parameters and r = 0.5 Rj up and 1 Rj up - We 
fit the predicted transit light curve to the combined pho- 
tometric data, and then compute the A% 2 between the 
constant flux fit and the predicted transit model. 

Our best-fit has a A\ 2 = —19.8, which we do not con- 
sider significant. We find similar or larger improvements 
in x 2 when we consider arbitrary phases for the tran- 
sit and when we consider "anti-transits" (signals with 
the same shap e as transits but cor responding to positive 
deviations, see lBurke et al. Il2006[ ). As before, these for- 
mally significant signals likely arise from systematics in 
the photometric data. We conclude there is no evidence 
for a transit signal in the combined data. 

Given that we do not detect a transit signature, we 
can also use this procedure to determine the confidence 
with which we can rule out transits of a companion with 
a given radius. This is just given by the fraction of the 
steps in the Markov Chain for which the companion tran- 
sits and produces a transit signature with a Ay 2 relative 
to the fixed constant flux greater than some threshold. 
We consider thresholds of A% 2 = 9, 16, and 25. The re- 
sulting cumulative probability distributions for a range of 
companion radii are shown in Fig. 111! Given the system- 
atics in the data, we consider thresholds of A\ 2 > 25 to 
be robust, and thus conclude that transiting companions 
with r > 0.75 Rj up are likel y ruled out at the ~ 9 5% con- 
fidence level. The models of Baraffe e t al. I (]2003t ) predict 
radii of > 0.8 Rj up for brown dwarfs of m p ~ 60 Mj up 
and ages of < 5 Gyr. Given the upper limit of the age of 
the primary of ~ 2 Gyr, we can therefore essentially rule 
out non-grazing transits of the companion. 

7. CONSTRAINTS ON SYSTEM GEOMETRY AND TIDAL 
ANALYSIS 

The rotational period of a star can be estimated from 
an empirical r elationship between the Rossby number 
and \ogR' HK (|Noves et al. 1 119841 ). The Rossby num- 
ber i?o = P TjT 1 , where P is the rotation period of 
the star and r c is the convective turnover time. In 
this work, we use the r elatio nship as quantified by 
lMamaiek~ fc Hillen brand I ([2008). We estimate the con- 
vective turnover time based on the rela t ionsh ip be- 
tween R and (B - V) from iNoves et all (|1984[ ). We 
find an expected rotational period based on the mea- 
sured log R' HK of P = 2.93 ± 0.37 days. The un- 
certainty in the period includes the uncertainty in the 
IMamajek fc Hillenbrand] (|2008D relationship, as well as 
an adopted uncertainty in R' HK of 0.2 x 10~ 5 based on 
the observed varia bility of the most active stars in the 
lLovis et al. I (|2011l ) sample. The latter uncertainty ac- 
counts for the fact that the IMamajek fc Hillenbrand 1 
(2008) relation applies to a time-averaged R' HK , while 
we have a single epoch measurement. The modest up- 
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per limit on the rotation rate of usin/ < 6 km s _1 , 
the fact that the estimated rotation rate from the 
iMamaiek fc Hillenbrand 1 (]2008l ) relation is close to the 
orbital period of the inner companion, and the expected 
equatorial rotation velocity of 22.68 ± 3.12 km s _1 if the 
primary was rotating at the inner companion's orbital 
period and had a stellar spin axis oriented edge-on, all 
suggest that the primary's spin axis is inclined relative 
to our line-of-sight, in which case the inclination is con- 
strained to be / = 15.0tg;2 degrees. 

We can use the reasonably strong upper limit on TYC 
2930 's photometric variability to place additional lim- 
its on the inclination of the star and its companion. 
The lack of photometric variability is somewhat surpris- 
ing, given the spectroscopic indications that the star is 
relatively active. An estimate of the expected photo- 
metric yariability can b e obta ined using the relation- 
ship in lHartman et al. I (|2009T ) between Rq and pho- 
tom etric amplitude. Based on the estimated Rq from 
the IMamaiek fc Hillenbrand! (|2008ft i? -log R' HK rela- 
tionship, we would ex pect a photomet ric amplitude of 
-0.8%. Figure 17 in lHartman et al. I (|2009t ) indicates 
that there exists substantial scatter about this relation, 
which is likely partially suppressed for values of Rq > 0.4 
due to incompleteness. Extrapolating the observed scat- 
ter at lower values of Ro (which are less affected by in- 
completeness) , we expect a photometric amplitude in the 
range of ~0.2-2%. Thus, the fact that the observed am- 
plitude is a least a factor of ~3 times lower (< 0.07%) 
suggests that this star is either surprisingly photometri- 
cally quiet given its spectroscopic activity indicators, or 
is it being viewed nearly pole-on. 

The expectation of a synchronous rotation rate for 
the primary is reinforced by consideration of the tidal 
evolution of the primary and secondary. To analyze 
this effect, we employed the "constant-phase-lag" (CPL) 
tidal model in which the location of the tidal bulge on 
the two bodies lies at a constant angle from the line 
connecting the centers of the two bodies. Frictional 
forces inside the two objects prevent perfect alignment, 
which leads to energy dissipation and transfer of angu- 
lar mome ntum. Hence the system will evolve with time 
(see e.g. iGoldreich fc Soterl 119661: iFerraz-Mello et al.~l 
2008). Two outcomes are possible: arrival at the "dou- 
ble synchronous" state, in which the obliquities are nor- 
mal to the orbital plane, and the two spin frequencies 
equal the orbital fr equency, or the two bodies merge 
(|Counselman II1973D. Here we use the CPL model pre- 
sented in IFerraz-Mello et al. I (120081) , with the numerical 
methods described in iBarnes et al. I (|2012t App. D). 

In the CPL model, the tidal effects scale with the 
"tidal quality factor" Q*. This parameter is poorly con- 
strained in stars an d brown dwarfs with values ranging 
from 10 5 -10 9 (e.g. ILin et al. I [19961: IMatsumura et al. I 
2008); therefore we explore the timescale to reach the 
double synchronous state in this range. The masses and 
radii of the two objects are also important, and although 
we have constraints on the two masses and the primary's 
radius, the secondary's radius is unknown. We adopt a 
radius of 1 Jupiter m ass, in line with theoretical models 
(jBaraffe et al. |[2003h . 

For compact binaries, tidal evolution during 
the pre-main sequence phase is also important 



(IZahn fc Bouchet I 119891: iKhaliullin fc Khaliullina I 
I2011HG omez Maqueo Chew et al. II2012D . since the radii 
are larger. As tidal effects scale with radius to the 
fifth power, the results can be dramatic. Most stellar 
binaries with periods less than 8 days are on circular 
orbits and probably nea r the double synchronous state 
IZahn fc Bouchetl (|l989l) . suggesting the TYC 2930 sys- 
tem could also have reached this state early in its history. 
On the other hand, the inner companion could have 
arrived at its orbit after the radial contractio n phase, 
perha ps via a gravitational scattering event (|Heggie I 
119751 ). or due to Koza i-Lidov interaction followed b y 
tidal circularization ([Fabrvckv fc Tremaine I 120071) . 
Therefore, in order to be as conservative as possible, we 
will ignore any evolution during the pre-Main Sequence. 
We therefore consider our timescales to reach double 
synchronization to be upper limits. 

In Fig. [12] the evolutions of the orbital period, stel- 
lar spin period and stellar obliquity, ij), are shown for 
three different cases of initial obliquities: ip = 0° (solid 
curves), tpo = 10° (dashed curves) and ipo = 45° (dotted 
curves). For each case, the initial spin period is 30 days, 
the secondary has a mass of 100 Mj up , the initial orbital 
period is 2.8 days, and the stellar tidal quality factor Q* 
is 10 . In each model, the orbital and spin periods reach 
~ 2.43 days and become locked after ~ 300 Myr. The 
obliquity can evolve slightly longer, but with negligible 
effect on the spins and orbit. Double synchronization re- 
quires ip = 0° and hence the timescale to reach that state 
is larger when the initial ip is nonzero. 

The initial spin period of the star is unknown, and 
therefore we tested values larger than 30 days. For the 
extreme case of a 120 days period and V'o = 45°, the time 
to reach double synchronization is about 3 times longer 
than the cases shown in Fig. [12] For this initial config- 
uration, an initial orbital period of 2.85 days produces a 
better match to the observed system. 

We expanded our analysis to a range of M 2 and Q* to 
determine which values predict the double synchronized 
state on a timescale less than or equal to the system age 
of < 2 Gyr. Initially, we set the primary spin period to 30 
days and ipo to 0° (as shown in Fig. [121 different choices 
do not affect our results) , the secondary's spin period was 
always tide-locked, and the initial orbit was circular with 
a period of 2.8 days. We then integrated forward a suite 
of configurations in the ranges 60 < Mi < 1000 Mj up 
and 10 6 < Q* < 10 9 and calculated the timescale to 
reach double synchronization. We considered this state 
achieved when the stellar spin period was within 10% of 
the orbital period. Fig.lTSIshows the timescale for double 
synchronization as a function of inner companion mass 
and Q*. For the minimum mass, double synchronization 
will occur within the nominal 2 Gyr lifetime if Q* < 
3 x 10 7 , while at M 2 = 300 M Jup , Q» must be < 4 x 10 8 . 
Our tidal analysis strongly favors the double synchronous 
state. 

8. A POSTERIORI DISTRIBUTIONS OF THE TRUE MASS 

The a posteriori distribution of the true mass of the 
companions given our measurements depends on our 
prior distribution for the mass of the companion, or, 
roughly equivalently, our prior on the mass ratio. As 
a rough illustration, if we assume a prior that is uniform 
in the logarithm of the true mass of the companion, then 
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the distribution of cos i will be uniform. Therefore, the 
median cosi ~ 0.5, and thus the median sini ~ 0.866. 
Thus we have, 

0.75 3/2 M? - Mj (Mi + Mjf = 0, (9) 

where j = 2,3 refers to the secondary and tertiary, re- 
spectively. For the secondary, A^2 = 1-6711 x 10 and 
Mi = 1.21 M ; the solution is M 2 ~ 0.0752 M Q or 
~ 78.75 M Jup . This result is roughly M 2 , m m/v / 07r5, 
but not exactly. For the tertiary, the median mass is 
M 3 - 0.48 M . 

More generally, for other priors, cos? is not uniformly 
distributed. We adopt priors of the form: 



where q is the mass ratio between the companion and the 
primary, and a = — 1 for the uniform logarithmic prior 
discussed above. To include this prior, we draw a value of 
cos i from a uniform distribution for each link the MCMC 
chain, but then weight the resulting values of the derived 
parameters for that link (i.e., the companion mass m) by 
q a+1 . In addition, based on the analysis of the relative 
photometry in §6.3[ we exclude transiting configurations 
of the secondary such that b = acosi2 P] -1 < 1. 

For q > 0, the a posteriori distribution does not con- 
verge, i.e., there is finite probability at infinitely large 
true masses. However, we can rule out nearly equal mass 
ratio main-sequence companions by the lack of any evi- 
dence of additional light beyond that from the primary, 
in particular from the lack of a second set of spectral lines 
and from the shape of the spectral energy distribution. 
As we describe below, both of these constraints place an 
upper limit on the mass of any main-sequence companion 
to < 1 M Q . We adopt this limit by giving zero weight 
to inclinations of the secondary and tertiary such that 
M2,3 > 1 M©. In doing so, we implicitly assume that 
neither of these companions are stellar remnants. 

The optical/near- infrared colors of the primary are 
well-fit by a model for the spectral energy distribution 
appropriate to a star with the effective temperature and 
surface gravity we measure from the spectrum (see Fig. 
[1]). A sufficiently massive main-sequence companion 
would contribute near-IR flux in excess of that predicted 
from the primary. We can therefore use the lack of ex- 
cess flux to constrain the mass of the secondary (and 
tert iary). We gene r ate a set of two-component SEDs us- 
ing iBaraffe et al~l (|1998l I2003D models, adopting stellar 
parameters for the primary from Table [3] and models cor- 
responding to tertiary masses from 0.2 to 1.2 M Q . We 
then fit the observed fluxes with these two-component 
SEDs, allowing distance and line-of-sight extinction to 
be free-parameters. For each tertiary mass, we calcu- 
late a x 2 = Xsed + xLu where xIeb is based on the 
observed fluxes and xlist i s based on the Hipparcos par- 
allax. We compute a A% 2 = xLomp ~ Xicomp. where 
Xicomp i s the x 2 f° r the single-component (i.e., primary- 
only) SED fit. We then use these A^ 2 to determine an 
upper limit of the tertiary mass by rejecting those ter- 
tiary masses that yield A^ 2 > 0, finding M3 < 1 M©. 
A similar mass constraint is derived by visual inspec- 
tion of the cross-correlation functions from the spectro- 



scopic observations. Assuming none of the components 
are evolved, mass ratios of ~ 0.8 can be excluded at the 
maximum RV separation between the primary and ter- 
tiary (~ 5 km s _1 ) due to a lack of asymmetry in the 
correlation peaks, once again constraining M3 < 1 M . 

Finally, as discussed in <|7l the rotational period of 
the primary P» as estimated from the R' HK index is 
within 1.5(7 of the period of the secondary, suggest- 
ing that the primary spin may be synchronized to the 
companion orbit. This hypothesis is corroborated by 
our tidal analysis, which suggests that for reasonable 
values of the stellar Q„ this system should reach a 
double-synchronized state in which the obliquities are 
aligned with the orbit within 100 Myr, which is con- 
siderably less than our estimate of the age of the sys- 
tem. Therefore, adopting this synchronized/coplanar 
assumption, we have I = i% and P* = P2. We can 
use these assumptions, combined with the constraints 
on wsin/ = 3.8+28 km s" 1 from the spectroscopic anal- 
ysis of the primary, to constrain i% and thus the true 
mass of the companion. For each link the MCMC 
chain, we use the value of i 2 , R\, and secondary pe- 
riod P2 to estimate vsml — vsmiz — {2%Rx/ P%) sin«2- 
We then multiply the weight of that chain by the ad- 
ditional factor exp[— 0.5(i>sin/ — 3.8 km s _1 ) 2 /cr^ sinJ ] 
where <J VS u\i = 2.8 km s _1 for i>sin/ < 3.8 km s _1 and 
CusinJ = 1.9 km s^ 1 for wsin/ > 3.8 km s _1 . 

Figs. Q31 [15l and [TS] show the resulting cumulative a 
posteriori distributions of the true masses of the com- 
panions under several sets of assumptions. These fig- 
ures illustrate the effect of adopting various priors and 
constraints on our inferences about the nature of these 
companions. Tables |6] and [7] list the median and 68% con- 
fidence intervals on the true masses of the companions. 
We will discuss each model in turn. 

Fig. [Ml shows the results for the secondary, for a = — 1 
(uniform logarithmic prior on q) and a = (uniform lin- 
ear prior on q), assuming the constraint from the lack 
of transits and the upper limit on the mass of M© from 
the lack and evidence of light from the companion. The 
constraint that M2 < 1 M Q makes little difference for 
a = — 1, because the companion is unlikely to be suffi- 
ciently massive to contribute a significant amount of flux. 
For a — 0, the prior on the companion mass is more heav- 
ily weighted to larger masses, and as a result this con- 
straint does affect the high-mass tail of the probability 
distribution; nevertheless, the inferred median masses in 
the two cases are similar. Under these assumptions, we 
would conclude that the secondary is most likely a low- 
mass stellar or brown dwarf companion seen at a mod- 
erate inclination of ~ 40 — 60° , with a median mass just 
above the hydrogen-burning limit M 2 ~ 0.08 — 0.1 M©. 

Assuming synchronization and coplanarity changes the 
conclusion about the nature of the secondary dramat- 
ically, given that the spectroscopically- measured v sin / 
suggests a relatively low inclination for the primary of 
~ 15°. Fig. [TBI shows the inferred cumulative distribu- 
tion for the secondary, for a — and a = — 1 and as- 
suming the upper limit on the mass. The results for the 
two priors are broadly similar: under the synchronization 
assumption the secondary is most likely a mid M dwarf 
with M 2 ~ 0.3 — 0.4 M Q with a nearly pole-on orbit with 
i 2 - 10-13°. 



10 



Scott W. Fleming et al. 



Fig. 1161 shows the cumulative distributions of the true 
mass of the tertiary. The mass of the tertiary is limited 
from below to be > 0.2 M by the measured mass func- 
tion, and from above to be < 1 M by the lack of light 
from the companion. The choice of prior on q has a rel- 
atively weak affect on the distribution of masses within 
this range. The tertiary is most likely to be an early M 
dwarf with a mass of ~ 0.5 — 0.6 M with a moderate 
inclination of is ~ 60 — 70°. 

Adopting the synchronization/coplanarity assumption, 
we infer that the secondary is seen nearly pole-on, 
whereas from the lack of evidence for additional light 
from a companion, we infer that the tertiary cannot be 
very massive and thus cannot have a low inclination. Fig. 
1171 shows a posteriori probability densities for the incli- 
nations of the tertiary and secondary under these as- 
sumptions. The 95% confidence level (c.l.) upper limit 
on the secondary inclination including the synchroniza- 
tion/coplanarity constraint and flux constraint is 23.6° 
(a = -1) and 21.2° (a = 0). On the other hand, the 
95% c.l. lower limit on tertiary inclination including the 
flux constraint is 30.1° for a = — 1 and 28.1° for a = 0. 
The difference in the orbital inclinations with respect to 
the sky plane \i\ — i 2 | is a lower limit to the true mutual 
inclination, and thus under these assumptions the orbits 
are misaligned at the > 95% c.l. 

The true mutual inclination A023 of the orbits is re- 
lated to the inclinations referenced to the sky plane by, 

cosA</>23 = cosi 2 cos«3 + sini 2 sini 3 cos(Af2 23 ), (11) 

where A51 2 3 = 51 2 — ^3 and flj is the longitude of the 
ascending node for companion j. Fig. 1181 shows the dis- 
tribution of mutual inclinations assuming the flux ratio 
constraint, a uniform distribution for A57i 2 , and the two 
different priors on a. These distributions are symmetric 
about A0i 2 = 90° because the observables (radial veloc- 
ity amplitude and potential transits) are invariant under 
the transformations ij — > ij + n. The 95% c.l. lower and 
upper limits on the true mutual inclinations are 31.1° 
and 180° - 31.1° = 148.9° for a = -1 and 28.8° and 
180° - 28.8° = 151.2° for a = 0. The probability that 
the mutual i nclination is greater than the Kozai angle 
(|Kozai I fl 9 6 2f) of 39.2° or less than the retrograde Kozai 
angle of 140.8° is ~ 89% (a = -1) and - 85% (a = 0). 

Fig. [19] shows the probability densities of the true mass 
of the secondary and tertiary, assuming the tidal syn- 
chronization/coplanarity and flux ratio constraints, for 
a = — 1 and a = 0. We infer that the primary and sec- 
ondary are both likely to be M dwarfs, with the tertiary 
likely to be somewhat more massive than the secondary. 

9. THE KOZAI-LIDOV MECHANISM APPLIED TO TYC 

2930 

One mechanism of forming short-period binary stars 
is the Kozai-Lidov mechanism, in which a hierarchical 
triple system exchanges angular momentum between the 
inner and outer orbits periodically if the mutual incli- 
nation between the inner and outer orbit s is 39.2° < 
Afe < 141.8° (|Kozai I Il96l ILidov I fl962b . These os- 
cillations drive e 2 to high and low values while decreas- 
ing and increasing A</> 23 , respectively. When combined 
with tidal friction, these oscillations can cause the semi- 
major axis of the inner companion (a 2 ) to decrease, as 
tidal friction removes energy and circularizes the or- 



bit (jMazeh fc Shaham I ["1979ft . In addition to short- 
period binaries in triple system s (Kisel eva et al. 1 119981: 
l Eggleton fc Kiseleva-Eggleton I 120011: iTokovinin et al. I 
2006; Fab rvckv fe Tremaine 1 120071 ) . this mechanism has 
been proposed as an explanation for the prevalence 
of "Hot J upiter" gas giant planets with periods of a 
few days dWu fc Murray 1 120031: iFabrvckv fe Tremaine I 
[20071 IWu et al. 1 120071 INaoz et al. 1 1201111 and the for- 
mation of blue straggle r stars in globular clusters 
(|Perets fc Fabrvckv 1 12001 . 

To constrain the initial parameters of TYC 2930 's 
progenitor system, we first assume that a 3 and e 3 re- 
main unchanged during the evolution of the system. We 
can then place an upper limit on the initial a 2 (a 2j i) 
using the stability criteria for triple systems given by 
(|Mardling fc Aarseth |[200l : 



03 

«2 



>Cf 



1 



M 3 



Ah+AhJ (1-ea) 



l + e 3 



(12) 



where C — 2.8 is an empirically-fit constant and / = 
1— — A</> 2 3 is an ad-hoc mutual inclination term. Given 
the most probable values for the parameters of TYC 
2930, and assuming that the initial A</> 2 3 ~ 90° such 
that the Kozai-Lidov mechanism is a significant effect, 
the initial 0,3/0,2 > 4.4, and thus the maximum stable 
02 is a 2 ,max = 0.96 AU. To replicate the observed sys- 



tem parameters, we then require 



peri 



< 3i?i = rtidc 



for efficient tidal circularization. Generally, the cos A0 23 
required for r per i < 3 Ri, cos A0 2 3 crit , is given by: 



cos A0 23c 



02.i 



6 r tidc 
5a 2 ; 



1/2 



(13) 

Thus if a 2 ,i = 0.3 AU then cos A0 23crit - 0.25. 

To simulate the Kozai-Lidov mechanism, we use a 
modified version of the A-body code FEWBODYE3 
(|Fregeau et al. Il200l . In Fig. [20l we show an example 
triple system which has initial parameters consistent with 
the constraints placed on the progenitor system of TYC 
2930: Mi = 1.21 M , M 2 = 0.34 M , M 3 = 0.48 M , 
a 2 ,i = 0.3 AU, a 3 = 4.27 AU, e 2 = 0.0, e 3 = 0.29, 
cj 2 = 0°, cj 3 = 0°, and cosA</> 23 = 0.1. The top panels 
present the evolution of a 2 , 03 and e^. The middle and 
bottom panels show the evolution of cos A(f>23 and 1 — e 2 , 
respectively. In the middle panel, the teal line designates 
cosA(/> 2 3 = 0. In the bottom panel, the red line shows 
where e 2 is high enough such that the radius of perias- 
tron of the inner binary (r por ;) is equal to 2 R . The 
example system goes through Kozai-Lidov cycles, bring- 
ing r pcr i < 2 R . If the effects of tides were included 
in our calculation, we would expect that a 2 and A0 2 3 
would decrease over many Kozai cycles. 

However, the commonly-used, quadrupole-order ex- 
pansion of the three-body Hamiltonian is insufficient 
to capture the secular dynamics of triple systems un- 
der the test particle approximation when e 3 is non-zero 
(iLithwick fc Naoz ll201ll : iKatz et al.l 120111: INaoz et al. I 
120111 ). The importance of the octupole-order terms 
relative to the quadrupole-order terms in the doubly- 



25 FEWBODY is 

http : 111 ewbody . sourcef orge .net/ 



available 



at 
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averaged three-body Hamiltonian is given by the param- 
eter 



Eoct — 



(M 


-MA 




\M 


+ mJ 





In the limit that (M 2 < Mi,M 3 ) and e 3 ^ 0, it is 
possible for the triple system to "flip", i.e., the sys- 
tem exhibits quasi-periodic cycles in A</> 23 through 0, 
and the tertiary passes between prograde an d retrograde 
(|Lithwick fe Naoz 1120111; iKatz et al. 1120111 ) . These flips 
occur even for small values of e oct (~ 10~ 3 ) in the test 
particle approximation, as long as the system is suffi- 
ciently inclined and the arguments of p eriastrons are cho- 
sen judiciously (see Figures 3 and 7 of lKatz et al. |[20TTI 
and iLithwick fc Naoz 11201 U respectively). 

Flips become increasingly common for triple sys- 
tems with larger values of e oc t, and they corre- 
spond to extremely large spikes in e 2 , such that 
(1 — e 2 ) ~ 10 -5 . Such qualitative behaviors (flips and 
eccentricity spikes) are referred to as the "eccentric 
Kozai-Lidov mechanism" , and can occur over a broad 
range of parameters wi t hin o ctupole-order calculations. 
IShappee fc Thompson I (|2012l ) have recently highlighted 
the mass dependence of the eccentric Kozai-Lidov mech- 
anism, and performed one of the first explorations of the 
eccentric Kozai-Lidov mechanism for triple stellar sys- 
tems. They found that many triple systems become 
tidally affected while on the main sequence as a result 
of this mechanism. The mechanism can even be effec- 
tive for triple systems that would not otherwise become 
tidally affected during normal Kozai-Lidov oscillations in 
quadrupole-order calculations. 

As an example of the eccentric Kozai-Lidov mecha- 
nism, we integrated the orbit of a system with M\ = 
1.21 M Q , M 2 = 0.34 M Q , M 3 = 0.48 M Q , a 2 = 0.6 AU, 
a 3 = 4.27 AU, e 2 = 0.5, e 3 = 0.35, u 2 = 0°, w 3 = 262°, 
and cos A023 =0.1. Fig. |2"T1 shows the resulting evolu- 
tion. There are many qualitative differences between the 
evolution of this system and that shown in Fig. [20l The 
most obvious of these are the periodic oscillations of e 3 , 
the long-term oscillations of the minimum cos A0 23 , and 
the nipping of the sign of cos A0 23 , which corresponds to 
large spikes in e 2 ~ 0.9995. 

What is the fate of systems that become tidally af- 
fected due to the eccentric Kozai-Lidov mechanism? Un- 
fortunately, no study has been performed to investi- 
gate the eccentric Kozai-Lidov mechanism with the in- 
clusion of tidal dissipation for triple stellar systems. It 
is possible that tidal friction will simply detune the ec- 
centric Kozai-Lidov mechanism, inhibiting the extreme 
spikes in eccentricity o bserv ed in our example system. 
However, iNaoz et al. I (|2011| ) investigated the eccentric 
Kozai-Lidov mechanism with tidal friction for a Jupiter- 
mass secondary in an octupole-order calculation. They 
demonstrate that the extremely high e 2 obtained during 
a flip can lead to large tidal effects, which rapidly cap- 
ture the planet in a short-period retrograde orbit. They 
coin this mechanism "Kozai capture" and claim it may 
explain the occurrence of retrograde, short-period, gas 
giant planets. A similar mechanism may operate for stel- 
lar triple systems, in which Kozai capture would cause a 
rapid dissipation of orbital energy from the inner binary 
during its evolution. This would result in a circularized 
inner binary with a semimajor axis of a few r per i, and 



thus provides a plausible scenario by which TYC 2930 's 
secondary star was driven to such a small semimajor axis. 

10. CONCLUSION 

We have discovered a short-period companion to TYC 
2930-00872-1 with a minimum mass below the hydrogen 
burning limit. Despite its relatively high transit prob- 
ability, we exclude any transits of the companion with 
high confidence using data from three ground-based tele- 
scopes. A long-term trend in the RVs indicate the pres- 
ence of a longer-period tertiary in the system. The ter- 
tiary's spectral lines are not detected in our spectroscopic 
data, its fluxes do not significantly contribute to our SED 
fitting, and direct imaging excludes stellar-mass, main- 
sequence companions out to projected separations of 30 
AU. Our spectra show the clear presence of Ca H and 
K core emission, but there is an unexpected lack of pho- 
tometric variability, and the measured v sin / is signif- 
icantly smaller than expected if the primary's rotation 
rate was tidally synchronized to the inner companion's 
orbital period. This suggests the primary's stellar spin 
axis is closely aligned to the line-of-sight. Given the age 
of the system, it is expected that the inner companion's 
orbital angular momentum vector is aligned with the stel- 
lar spin axis, therefore its line-of-sight orbital inclination 
is low, and its true mass is likely to be stellar. The ab- 
sence of any detected signal from either component in the 
spectra and SED place an upper mass limit of ~ 1.0 M©, 
if they are not stellar remnants. Assuming the tertiary 
is not a remnant, the upper mass limit places a lower 
limit to its line-of-sight inclination, which results in a 
significant mutual inclination between the secondary and 
tertiary. Such mutual inclinations are expected if the sys- 
tem's dynamical history was driven by the Kozai-Lidov 
mechanism. Long-term RV monitoring of the outer com- 
panion to obtain reliable orbital parameters will greatly 
improve the constraints that can be placed on the mu- 
tual inclination between the secondary and tertiary. Fur- 
thermore, high signal-to-noise ratio spectroscopic obser- 
vations that could detect the presence of the (presum- 
ably) M dwarf companions would also allow for masses 
to be assigned to both objects and improve the inclina- 
tion constraints. 
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Table 1 

MARVELS RV - TYC 2930 



HJDutc 


RV 


QF-Scaled ctrv 




(km s- 1 ) 


(km s- 1 ) 


24o4o4o . 8oy4u 


2.571 


0.129 


2454844. 8259o 


-8.036 


0.099 


2454845. ooooo 


-6.326 


0.095 


2454840.82454 


1.039 


0.077 


2454847.77335 


-13.549 


0.083 


2454800.71528 


-7.996 


0.087 


o /i cr /i C7/i 7J007 

2454874. 74o27 


-10.879 


0.118 


2454876.77555 


-12.778 


0.128 


2455135.87308 


9.021 


0.120 


24551oo. 84418 


-7.602 


0.105 


2455137.85985 


7.069 


0.165 


2455138.88972 


-1.846 


0.120 


2455139. 78oo8 


-3.523 


0.110 


24551oy.yoyyi 


0.506 


0.099 


245514a. 8o472 


-4.111 


0.102 


2455144. 8o2o2 


1.319 


0.084 


O/icrcri/icr oca 7n 

2455145. 8o47U 


5.141 


0.107 


2455171.89913 


7.588 


0.094 


a et cr 1 to oc7i 

2455172.85712 


-0.125 


0.118 


2455200.91808 


5.337 


0.077 


2455254.77584 


10.438 


0.085 


2455258. /51 (4 


-4. / 26 


0.104 


2455281.65360 


10.108 


0.141 


2455466.86921 


1.955 


0.137 


2455466.91135 


0.996 


0.105 


2455487.84516 


7.664 


0.077 


2455488.86530 


-0.942 


0.101 


2455489.87545 


-0.646 


0.095 


2455494.90526 


3.061 


0.099 


2455500.90889 


1.409 


0.092 


2455521.82966 


6.852 


0.107 


2455522.85494 


-0.199 


0.117 


2455550.88590 


4.448 


0.195 



Table 2 

SARG RV - TYC 2930 



HJDutc RV o- RV 

(km s _1 ) (km s _1 ) 



2455436. 


71196 


2 


.568 


0. 


.026 


2455460. 


73287 


-3 


.624 


0. 


015 


2455460. 


74384 


-3 


.400 


0. 


018 


2455460. 


75505 


-3 


.196 


0. 


021 


2455495. 


57813 


5 


.964 


0. 


.034 


2455495. 


61314 


5 


.550 





.020 


2455495. 


69113 


4 


.642 


0. 


.024 


2455495. 


71406 


4 


,247 


0. 


.023 


2455516. 


.55620 


-5 


.453 


0. 


.022 


2455516. 


63469 


-3 


.758 


0. 


012 


2455553. 


50056 


4 


,252 


0. 


.014 


2455553. 


66834 


5 


.920 





013 


2455580. 


.42134 


5 


.869 


0. 


.012 


2455580. 


46981 


6 


.010 


0. 


.012 


2455580. 


59158 


5 


.770 





.014 


2455666. 


.42077 


-8 


.355 





014 


2455698. 


35676 


-12 


.133 


0. 


016 


2455791. 


70794 


3 


,088 





.012 


2455844. 


61269 


-8 


.335 


0. 


015 


2455844. 


74543 


-5 


.583 


0. 


012 
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Table 3 

Stellar Properties - TYC 2930 





Parameter 




Value 




±1<7 


fy T o nnn [ dpf t a 

^.12000 


93.880921 


0.000004 


"J2000 


+39.931826 


0.000005 


FUV b 


19.815 


0.195 


NUV b 


14.34 


0.01 


B (HAO) 


10.365 


0.023 


V (HAO) 


9.842 


0.018 


J (2MASS) 


8.770 


0.029 


H (2MASS) 


8.539 


0.047 


K s (2MASS) 


8.458 


0.023 


WISE 3.4 //m 


8.380 


0.024 




o.oyz 


u.uzo 


WISE 12 pim 


8.329 


0.029 


WISE 22 pirn 


8.201 


0.245 


(mas yr _1 ) c 


3.42 


2.05 


/ — l\c 

fi s (mas yr J 


a r fin 


1 1 Q 

1. lo 


Parallax n (mas) c 


/ .ID 


1.01 


A v (SED) 


0.33 


0.06 


T cff (K) 


6427 


33 


log(s [cm s *]) 


4.52 


0.14 


[Fe/H] 


-0.04 


0.05 


£t (km s x ) 


1.40 


0.05 


v sin/ (km s _1 ) 


3.8 


+ 1.9 

-2.8 


M. (M@) 


1.21 


0.08 


R* (R ) 


1.09 


+0.15 
-0.13 


RPMj 


2.14 





a Tycho-2 Catalog l lH0g et al.~ll2000l) 
b GALEX IMartin et al. ||2005|) 
c Ivan Leeuwen I 1120071) 



Table 4 

IAC and BPG Stellar Parameters 





IAC 


BPG 


Parameter 


SARG 


ARCES 


SARG 


ARCES 


Fe I lines used 


173 


172 


60 


67 


Fe II lines used 


21 


25 


8 


9 


Toff (K) 


6456 ± 49 


6413 ± 41 


6406 ± 110 


6415 ± 76 


log (g [cm s !]) 


4.68 ± 0.27 


4.53 ±0.21 


4.47 ± 0.26 


4.44 ±0.21 


[Fe/H] 


-0.02 ±0.09 


-0.01 ±0.07 


-0.13 ±0.10 


-0.03 ±0.07 


ft (kms" 1 ) 


1.296 ± 0.076 


1.464 ± 0.059 


1.44 ±0.20 


1.34 ±0.12 






Table 5 








Orbital Parameter 


s For The Inner Companion 






Parameter 


Median 


<T (Low) (T 


(High) 




P (days) 2.430420 


0.000006 0.000006 




Ktkms- 1 ) 8.723 


0.009 


0.009 


T c (HJD UTC - 2450000.0) 4842.2640 


0.00187 0.00187 




e 0.0066 


0.0010 


0.0010 




141.948 co 142 


7 


7 




7o,inst (km s 1 " 


1 -6.642 


2.106 


0.904 




7off (km s- 1 


1 -2.890 


2.102 


0.900 




70 (km s — 1 / 


1 35.751 


0.285 


0.285 
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Table 6 

A priori secondary mass. 



a a key b Median M 2 , 68% Median i 2 , 68% Transit Prob. 

-1 1 0.080^°°^ 56.0i^-j 0.133 

1 0.098t°;J6° 42.4127-2 .086 

-1 2 0.338±°; iJq 12.6±4 j < 0.001 

2 0.416to^5° 10.6i|;g < °- 001 



a Prior of the form: N/dq oc q a 

^ 1— Flux Ratio & Transit Constraint. 2— Flux Ratio, Transit, and 
Synchronization / Coplanarity. 



Table 7 

A priori tertiary mass. 



a Median M 3 , 68% Median z 3 , 68% Transit Prob. 

-1 0.480l°;^ 65.4i^;2 0.001 

0.553l°;^3 61.5±2§"_4 0.001 



a Includes Flux Ratio Constraint 
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Figure 1. NextGen model (solid line) compared to the observed broadband fluxes of the host star. Blue points represent the expected 
fluxes in each band based on the model, red horizontal bars are the approximate bandpass widths, and red vertical bars are the flux 
uncertainties. The T e g, log (g) and [Fe/H] are fixed at the spectroscopically-determined values, while Ay is allowed to float. No evidence 
of IR excess is detected, while there is potentially some GALEX FUV excess indicating elevated levels of stellar activity. 



18 



Scott W. Fleming et al. 




i i i i i i i i i i i i i i i i i i i i i 
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T ef f 

Figure 2. HR diagram based on Yonsei-Yale stellar evolution models (Dcmarque ct al. 2004). The solid track is for the best-fit stellar 
parameters, while the two dashed tracks represent the \-a uncertainties. The blue dots represent star ages in Gyr. TYC 2930 (red point) 
appears to lie on the main sequence. 
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Figure 3. Ca II K line of TYC 2930 compared to standard stars. TYC 2930 has clear core emission indicating that the host star is active 
compared to other stars with similar stellar parameters. 
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Figure 4. Phase-folded MARVELS (blue) and SARG (red) RVs of MARVELS-2b. The orbit of the long-period companion and the 
systemic velocity (7o,inst = —6.642 km s _1 ) have been removed in the top panel. The bottom panel shows the residual RVs after removing 
the short-period companion's orbit and systemic velocity. The RV uncertainties are not visible at this scale. 
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Figure 5. Unfolded MARVELS (blue) and SARG (red) RVs after removing the short-period companion's orbit and the systemic velocity 
(70,inst = —6.642 km s _1 ). The best-fit model of the outer companion is overplotted as the black line. The residuals of the combined, 
two-companion model are shown in the bottom panel. 
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Figure 6. Lucky images of TYC 2930 using the best 15 and 85% of the frames for the 2010 and 2011 observations. The CCD gain, which 
changed during the 2011 observations, is labeled by g in the images. No tertiary companion is detected. 
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E 
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Figure 7. Keck AO image of TYC 2930. No stellar companions are detected at the 3-<r level beyond ~200 mas ~ 30 AU, corresponding 
to a brightness limit of Am ~6. 



24 



Scott W. Fleming et al. 




20 30 40 50 60 70 80 



projected separation (AU) 



Figure 8. Dctcctabi lity (contrast cu rve) for the Lucky Imaging and Keck AO images of TYC 2930. Contrast levels are converted to 
masses based on Baraffc ct al. (2003) models for the Keck band and [Girardi ct al. (2002) isochrones for the Lucky Imaging band. A 
separation of 50 AU is ~ 350 mas at TYC 2930's Hipparcos-bascd distance of ~ 140 pc. 



MARVELS-2B 



25 




4000 4500 5000 5500 

HJD-2450000. 
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1 10 100 1000 10 4 

Period (days) 

Figure 9. Top panel: Cleaned relative photometry of TYC 2930 from WASP (blue), KELT (red), and Allegheny (green). Bottom panel: 
Lomb-Scargle periodogram of the combined photometric data. While a large number of strong peaks are visible, we do not regard these as 
significant. There is no strong peak at the period of the secondary (vertical red dotted line) or within the estimated 2 a period range of the 
primary's rotational period (vertical blue dotted lines). The inset shows detail of the most significant peak in the combined periodogram 
(grey), as well as the periodograms for just the WASP (blue dotted), KELT (red short dashed), and Allegheny (green long dashed) datascts. 
The peak in the combined dataset arises almost exclusively from the WASP data, and is not confirmed by the KELT data. 
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Figure 10. Relative photometry folded at the period of secondary, and binned 0.05 in phase. Phase zero corresponds to the expected 
time of conjunction (and so of transits for the appropriate inclinations). Black points are the combined data, blue are WASP, red are 
KELT, and green are Allegheny. The grey curves show the expected transit signatures for a companion with radius of 0.5 Rj up (dashed) 
and 1 Rj up (dotted), assuming an edge-on inclination and the median estimated values of the primary mass and radius. The solid curve 
shows the expected signature of ellipsoidal variability assuming an edge-on companion. 
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Figure 11. Probability that transits of a companion are excluded at levels of Ax 2 = 9, 16,25 based on the analysis of the combined 
WASP, KELT, and Allegheny photometric data sets, as a function of the radius of the companion. Transits of companions with radius 
r > 0.7 Rjup can be excluded at the 95% confidence level. 
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Figure 12. Tidal evolution of the TYC 2930 system. Curves represent assumed initial stellar obliquity values of ipQ = 0° (solid), 10° 
(dashed), and 45° (dotted). Top: Evolution of the orbital period. Middle: Evolution of the stellar spin period. Bottom: Evolution of the 
stellar obliquity. 




Figure 13. Time (years) for TYC 2930 to reach double synchronization. The star initially has no obliquity and a rotation period of 30 
days. The companion is tidally locked. No non-tidal effects, such as an early epoch of radial contraction, are included. 
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Figure 14. The cumulative a posteriori probability of the true mass of the secondary, for two different priors on the mass ratio q: 
dN/dq oc q^ 1 (uniform in logg), and dN/dq = constant (uniform in q). Companion masses > 1 Mq and transiting configurations are 
excluded. 



MARVELS-2B 



31 



x 
cd 
x 
o 

S-H 

> 



3 



Secondary, Synchronized and Coplanar 



0.8 



0.6 



0.4 



0.2 - 



i 1 1 1 — i — i r 



t 1 1 1 — i — i — r 



dN/dq «q~ 



dN/dq = constant 




0.01 




0.1 

Mass (M G ) 



_l I I I L 



100 

Mass (M Jup ) 



1000 



Figure 15. The cumulative a posteriori probability of the true mass of the secondary under the same set of assumptions as Figure H"4l 
but assuming that the primary and secondary are synchronized and the stellar obliquity has been damped to zero, and thus the inclination 
of the star / as determined from the observed velocity broadening sin I and stellar rotation period P„ is equal to the inclination of the 
orbit of the secondary i. 
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Figure 16. The cumulative a posteriori probability of the true mass of the tertiary under the same set of assumptions as Figure H"4l except 
that transiting configurations are not excluded. 
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Figure 17. The a posteriori probability densities for the inclinations of the tertiary and secondary. These curves have been arbitrarily 
normalized by their peak probability density. The black lines show the probabilities for the secondary assuming the tidal synchroniza- 
tion/coplanarity, for two different priors on the mass ratio q, dN/dq oc q (uniform in logg, solid), and dN/dq = constant (uniform in q. 
dotted). The grey lines show the probabilities for the tertiary, for the same priors on the mass ratio. In all cases, we exclude companion 
masses > 1 Mq. 
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Figure 18. The a posteriori probability density (black lines) and cumulative probability (grey lines) for the true relative inclination of 
the secondary and tertiary A</>23. The probability densities have been arbitrarily normalized by their peak probability density. The solid 
lines show the probabilities assuming a prior on the mass ratio q of dN/dq oc q^ 1 (uniform in logq), where as the dotted lines show dN/dq 
= constant (uniform in q). In both cases, tidal synchronization/coplanarity for the secondary has been assumed, and companion masses 
> 1 Mq have been excluded. The curves are symmetric about A0 = 90° because translating the inclinations of the secondary and tertiary 
by 7r results in the same observables. 
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Figure 19. The a posteriori probability density for the true masses of the tertiary and secondary. These have been arbitrarily normalized 
by their peak probability density. The black lines show the probabilities for the secondary assuming the tidal synchronization/coplanarity, 
for two different priors on the mass ratio q, dN/dq oc q^ 1 (uniform in log?, solid), and dN/dq = constant (uniform in q, dotted). The 
grey lines show the probabilities for the tertiary, for the same priors on the mass ratio. In all cases, companion masses > 1 Mq have been 
excluded. 
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Figure 20. Example of the normal Kozai-Lidov mechanism using parameters consistent with those of TYC 2930's progenitor system. 
Top: Evolution of a2 (dotted black), 03 (red dotted), and e$ (blue dotted). Middle: Time evolution of cos A<p23- The teal line designates 
cos A023 = 0. Bottom: Time evolution of e2, the red line denotes r por i < 2 Kq. Integrations after the first eccentricity maximum are for 
illustration purposes only, since we do not include the effects of tides here. 
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Figure 21. Similar to Fig. 1201 but for the eccentric Kozai-Lidov mechanism. Note the qualitative differences between the two cases: the 
changing sign of cos A<j>23 and the extreme eccentricity spikes. 
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